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The leading and next-to-leading threshold logarithms of the QCD corrections to 
electroproduction of heavy quarks in single-particle inclusive kinematics are re- 
summed to all orders in perturbation theory. The resummed cross-section is used 
to derive the NLO and NNLO results near threshold and their numerical impact 
on the charm structure function is studied. 



The deep-inelastic production of heavy quarks is an important reaction, 
which is most prominently used for the extraction of the gluon density in the 
proton. Experimental data on charm production in lepton proton collisions 
have been published by e.g. EMC 111, HI and ZEUSB, while calculations up to 
next^to-leading order (NLO) provide solid foundations for a theoretical descrip- 
tionu. As a matter of fact, the charm structure function is quite sensitive to 
the effects of partonic processes close to the charm quark pair threshold. This 
is due to the large gluon density g{x, ^) at small momentum fractions x, which 
enhances the contribution of soft gluon emission near the elastic limit □. In 
this kinematical region, the QCD corrections are dominated by large Sudakov 
double logarithms, which have to be resummed to all orders of perturbation 
theory. This task has been performed up to next-to-leading logarithmic (NLL) 
accuracy □. The technology and some results are briefly presentend here. 

We study electron proton scattering with the exchange of a single virtual 
photon, = —q^, and a detected heavy quark in the final state, thus 

j{q) + P{p) q{pi) + X , (1) 

where X denotes any additional hadrons in the final state and = w? . The 
Mandelstam invariants, S' — [ji ^ q)^ + Q"^ ,Ti = {p — piY — m? and Ui = 
{q ~ PiY ~ can be used define Si = S' + Ti + Ui, which vanishes near 
the partonic threshold. The double differential heavy quark structure function 
dF2 associated to the process eq.(||) may be written as 



d^F2{x,Si,Ti,Ui,Q) 



X! J ^ <l>i/p{y,^J')^2^(^^,S4,tl,Ul,Q, 



M , (2) 



dTidUl 

'=q.q,s ax 

where a = 1 + Am? /Q^ . The 4>i/p denote parton distributions in the proton 
at mometum fraction y and MS-mass factorization scale ji. The functions uj2i 
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describe the underlying hard parton scattering processes and depend on asso- 
ciated partonic Mandelstam variables s', ii, ui and S4, which are derived from 
eq.(|^) after replacing the proton P by a parton of momentum k = {x/y)p. At 
n-th order in perturbation theory, the gluonic hard part Lj2g in eq.(H) typi- 
cally depends on singular distributions a"[ln^"~^(s4/m^)/s4]_(-, that have to 
be resummed. Light initial state quarks will be neglected. 

In order to perform this resummation of threshold logarithms, the phase 
space near the elastic edge is decomposed into various regions, each depend- 
ing on an infrared safe kinematical weight ~ 0. For dF2 in eq.(H), this 
decomposition implies S^jiri} = wi + lUg, with w\ = (\ — x j y)[~u\) j rr? and 
Ws — S/^/m? . Then, general argumentsl3 allow for a refactorization of dF2^ 



dTidUi 



H2g{S,Ti,Ui) J dwidwsS (^-^ ~ 



Wl 



(3) 



separating the partonic degrees of freedom for Si/m? into different func- 
tions of individual weight w. The hard function H2g summarizes off-shell short- 
distances corrections, while ipg/g accounts for effects of collinear gluons and the 
soft function S for soft, long- wavelength gluons. Replacing the proton in eq.(||) 
by a gluon and taking Laplace moments, f{N) — dw exp[—Nw]f{'w), gives 



02g{N,h,Ui) = H2g{S,Ti,Ui) 



^g/giN,p-C) 

<i>g/g{N,fl) 



WCm) +0{1/N). (4) 



where (j^g/g is the usual MS-distribution from mass factorization. In moment 
space, the Sudakov logarithms appear as factors a"ln^"^W, with i — 0,1 for 
NLL accuracy, and the A''-dependence in eq. exponentiates for each function 
individually. All leading logarithms (LL) are exclusively contained in ipg/g, 
which is^a gluon distribution at fixed energy defined as an operator matrix 
elemental. It depends on a time-like vector (, ("^ — 1, that fixes the kinematics 
by projecting on the proper energy fraction p • of the proton. At order ag in 
D — 4, — 2e dimensions and axial gauge n ■ A ^ 0, n ~ ijjg/g is given by 



i'g/g{w,p ■ C = —Ca\— - + + 



U7 + \ /i^ 



(5) 



where the remaining collinear poles are cancelled by (j)g/g. The LL logarithms 
in i^g/g are resummed in analogy to the Drell-Yan process i, while all scale 
dependence of ipg/g ^-nd (t>g/g is governed by renormalization group equations 
(RGE) with anomalous dimensions 7^ and 7g/gl3. 
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The soft function S requires renormalization, since it is defined as a com- 
posite operator, that connects Wilson Hnes in the direction of the scattering 
partonsIM Its RGE, fj.{d /dfi) In S{N) = -2Rers, resums all logarithms in S 
and its gauge dependence cancels precisely the corresponding terms in V-'g/g- 
The soft anomalous dimension Ts is a one-dimensional matrix in colour space 
and to order calculated to be 



with f3 = ^l-Am^/s and = (1 - 2 mVs){ln(l - /3)/(l + /3) Itt}//?. 

The final result for the hard scattering function uj2g in moment space 
resums ajiJarge logarithms in single-particle inclusive kinematics up to NLL 
accuracy E3. Combining the resummed ipg/g with the integrated RGE for S, 
we obtain for (I'2g 

Q2g{N,h,ui)=H2giS,Ti,Ui) (7) 





p-C/jv 



1 



-Nw\ 



X 



exp{2 1 ^Rer5(a.(A2))}exp{-2y'^(7g(a,(A2))-7g/g(«.(A2)))}. 



The first exponent gives the leading A^-dependence of the ratio 'ipg/g/4'g/g 
with V(g){as) = 2CAas/Tr, A^g){as) = CAias/T^) + {CAK/2){as/Tt)^ and 
K = (67/18 - - 5/9nf, the latter ones being the well- known lid. 

Using the resummed result for L02g in eq.(H) as a generating functional, 
we reexpand 0)29 to NLO and NNLO. Inverting the Laplace transform and 
integrating over the phase space in cq.(|2|) leaves us for F2 with 

F2ix,Q)^^^^fl^ j d2/5(y,A^)E(47ra.(M^)f ^4';')(,7,01n'4-(8) 

where Cj^''"' are the standard gluon coefficient functionsi, = (^/4)(y/a; — 1) — 1 

and ^ = Q"^ jrr?. In fig.||on the left, we plot Cj^'^^ as a function of 77 for ^ = 4.4 
and compare the LL and NLL approximations up to NNLO. On the right, 
we plot the charm structure function F| for x = 0.1(0.01), = lOGeV^ as a 
function of a yma^-cut on the integral in eq.(@) to estimate threshold sensitivity. 
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Figure 1: Left: The jy-dependence of the gluon coefficient functions; exact results for 
C2j'°' and Cjg"'' (solid lines); approximate LL and NLL results for (cjg'"') (dotted 
and dashed lines). Right: Threshold dependence of the charm structure function F2 
at NLO with the CTEQ4M gluon distribution and = m = 1.6GeV; exact results 
(solid lines); approximate LL and NLL results (dotted and dashed lines). 



In general, at order as the NLL logarithms approximate very well, for 
moderate values of Q^, better than the LL logarithms, while the result for 
the gluon coefficient function S^^^^ represents the best present estimate of the 
NNLO corrections. An extensive derivation of the result eq.(||) for Cjig andj-a 
detailed study of its phenomenological consequences will be given elsewhereB. 
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